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The quantum partition function at finite temperature requires computing tlie trace of the imagi- 
nary time propagator. For immerical and Monte Carlo calculations, the propagator is usually split 
into its kinetic and potential parts. A higher order splitting will result in a higher order convergent 
algorithm. At imaginary time, the kinetic energy propagator is usually the diffusion Greens function. 
Since diffusion cannot be simulated backward in time, the splitting must maintain the positivity of 
all intermediate time steps. However, since the trace is invariant under similarity transformations 
of the propagator, one can use this freedom to "correct" the split propagator to higher order. This 
use of similarity transforms classically give rises to symplectic corrector algorithms. The split prop- 
agator is the symplectic kernel and the similarity transformation is the corrector. This work proves 
a generalization of the Sheng-Suzuki theorem: no positive time step propagators with only kinetic 
and potential operators can be corrected beyond second order. Second order forward propagators 
can have fourth order traces only with the inclusion of an additional commutator. We give detailed 
derivations of four forward correctable second order propagators and their minimal correctors. 



I. INTRODUCTION 



The quantum partition function requires computing the trace 

Z = Tr{p) = Tr{e-^"), (1.1) 

where p is imaginary time propagator, [3 — 1/ (ksT) is the inverse temperature and H = T+V is the usual Hamihonian 
operator. Although specific forms of the kinetic and potential energy operators will not be used in the following, it 
is useful to keep in mind the many-body case where T = (— ?i^/2m) X]i=i ^ ~ ^i<j ^(^u )- numerical or 

Monte Carlo calculations, the imaginary time propagator is first discretized as 



^-(3{T+V) ^ 



(1.2) 



where e = ~Aj3 = —[3/n, and the short-time propagator e^^"^'^^' is then approximated in various ways. One systematic 
method is to decompose, or split the short-time propagator into the product form 



e 



N 

e(T+y)^-Qgt.sTg..sy^ (1.3) 



with coefficients {ti,Vi\ determined by the required order of accuracy. For quantum statistical calculations, since 
(r'l e*'*^-^! r) cx e~^'' """^ /(2ti^/3) jg diffusion kernel, the coefficient U must be positive in order for it to be simulated 
or integrated. If ti were negative, the kernel is unbounded and unnormalizable, and no probabilistic based (Monte 
Carlo) simulation is possible. However, as first proved by Shengi, and later by Suzuki—, beyond second order, any 
factorization of the form H1.3|) must contain some negative coefficients in the set {ti, Vi\. Goldman and Kaper^ further 
proved that any factorization of the form H1.3|l must contain at least one negative coefficient for both operators. Thus, 
despite myriad of factorization schemes of the form (|1.3(l proposed in the classical symplectic integrator literature^iSiSii, 
none can be used for doing quantum statistical calculations beyond second order. It is only recently that fourth order, 
all positive-coefficient factorization schemes have been found^-^ and applied to time-irreversible problems containing 
the diffusion kerneliSiiiiiSiiiii. In order to bypass the Sheng-Suzuki's theorem, one must include other operators in 
the factorization 1)1.31) . such as the double commutator [V, [T, F]], where [A, B] = AB — BA. 

In computing the quantum partition function Z, only the trace oi p = e~^^ is required. Since the trace is invariant 
under the similarity transformation 

p = SpS~\ (1.4) 

one is free to use any such p to compute Z. This is immaterial if p is known exactly. However, if the short-time 
propagator is only known approximately, then one may use a clever choice of S to further improve the approximation. 
This is a well known idea in many areas of physics. For example, to calculate the exact quantum many-body ground 
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state using the Diffusion Monte Carlo algoritfim, one can choose S = (j)o, where (po is a known trial function close 
to the exact ground state. This is the idea of "importance sampling" as introduced by Kalos et alm^. Its operator 
formulation as described above has been implemented by Chinifi some time ago. Similar ideas have been used to 
improve path- integrals, as detailed by Kleinert^^. If the short time propagator is approximated by the product form 
(|1.3() . the error terms can be calculated explicitly and eliminated by S. When implemented classically, these are 
known as symplectic "corrector" , or "process" algorithmsiSiiSiSfliSiiSSiSi. In this context the propagator p is the kernel 
algorithm and S is the corrector. Since S disappears in the calculation of Z , there is no restriction on the form of S. 
If S were also expanded in the product form H1.3|) . there is no restriction on the sign of its coefhcients. This suggests 
that there may exist a product form H1.3() of p with only positive coefficients such that its trace is correct to higher 
order. This would not be precluded by the existing Sheng-Suzuki theorem. 

In this work, we show that this is not possible. If p is approximated by the product form with positive 

coefficients {ti}, then p cannot be corrected by S to higher than second order. The proof of this generalizes the 
Sheng-Suzuki theorem. The corrected propagator p can be fourth order only if additional operators, such as [V, [T, V]], 
are used in the splitting of p. By understanding the "correctability" requirement, we can systematically deduce the 
four fundamental correctable second order propagators and their correctors. 

In the following Section, we recall some basic results of similarity transforms. Beyond second order, only a special 
class of approximate p satisfying the "correctability" condition can be corrected to higher order. In Section III, we 
compute the explicit form of the error coefficients required by the correctability criterion. In Section IV, we show that 
this requirement cannot be satisfied for propagators of the product form with only positive {ti} coefficients. In 
Section V, based on our understanding of the correctability restriction, we deduce all four second order correctable 
propagators and their minimal correctors. Some conclusions are given in Section VI. 



II. SIMILARITY TRANSFORMS AND THE CORRECTABILITY CRITERION 



Similarity transforms on approximate propagators of the product form (|1.3() have been studied extensively in the 
context of symplectic correctorsiSiiSiSfliSiiS^. However, not all use the language of operators and some are specific to 
celestial mechanics. Here, we recall some elementary results and establish the fundamental correctability requirement 
in the context of quantum statistical physics. 

Since 

= [5e^(^+^)5-i]", (2.1) 

it is sufficient to study the similarity transformation of the approximate short-time propagator pA- Let pA approxi- 
mates e^'^+^^ in the product form such that 

p^ = ne*-^^e'"-^=e^^-, (2.2) 

where Ha is the approximate Hamiltonian 

Ha=T + V + eierv [T, V])+e^{ cttv [T, [T, V]] + evrv [V, [T, V]]) + 0{e^) (2.3) 

with error coefficients ctv, &ttv , &vtv determined by factorization coefficients {ti, Vi{. The transformed propagator 
is 

PA = SpaS-' = Se'"-S-' = e'iSH.'^-') ^ ^eS.^ ^2.4) 
where the last equality defines the transformed approximate Hamiltonian Ha- If now we take 

S = exp[£C] (2.5) 
where C is the to-be-determined corrector, then we have the fundamental result 

Ha = e^^ff^e-^c^ ^Ha+ e[C, Ha] + ^e^[C, [C, Ha]] + ^e^[C, [C, [C, Ha]]] + ■■■. (2.6) 

Let's first consider the case where the product form 1)2. 2|l for Ha is left-right symmetric, i.e., either ti — and 
Vi = VN^i+i, — tjsi-i+i, or Wat = and Vi = v^^i, ti = iTv-i+i. In this cases, the propagator is reversible, 
Pa{£)pa{^£) = 1, and Ha{s) is an even function of e with ctv = 0. In this case 

Ha - Ha+s[C,Ha] + ---, 

= T + V + e^{eTTv[T,[T,V]]+evTv[V,[T,V]])+s[C,T + ¥] + ■■■, (2.7) 



and one immediately sees that the choice C = eCi with Ci = ctv [T, V] would eliminate either second order error 
term with ctv = ^ttv or ctv = ^vtv- So, if Hj^ is constructed such that 

gttv — ^VTV (2-8) 

then both can be simultaneously eliminated by the corrector. This is the fundamental "correctability" requirement for 
correcting a second order pA to fourth order. This observation can be generalized to higher order. At higher orders, 
Ha will have error terms of the form [T, Qi] and [V, Qi] where Qi are some higher order commutator generated by T 
and V. If Ha is of order 2n in e, then Ha can be of order 2n + 2 only if Ha^s error coefficients for [T, Qi] and [V, Qi] 
are equal for all Qi's. This fundamental corrector insight is often obscured by the more general case where odd order 
errors are allowed. 

Sheng^ and Suzuki^ independently proved that no pA of the form l|2.2ll can have positive coefficients ti beyond 
second order. More precisely, if pA is of the product form H2.2(l with positive ti's such that ctv — 0, then cttv and 
evTV cannot both be zero. We will prove a more general theorem that the product form (|2.2(l with positive tiS such 
that Ctv = cannot be corrected beyond second order, i.e., ettv can never equal to evTV- From this perspective, 
the Sheng-Suzuki theorem is a special case where the common value for both coefficients is zero. 

In the general case where eyy ^ 0, we have 

Ha =T + V + e{eTv [T, V])+e\ cttv [T. [T, V]] + evrv [V, [T, V]] ) 

T + V] + e^TvlC, [T, V]] + \e\C, [C, T + V]]+ 0{e^). (2.9) 

Since \ctT + cyV., T ^V\~ (ct — cv)[T, V], the liner term in e can be eliminated if we choose C — Cq = ctT + cyV 
such that 

{cT'Cv)^-eTv (2.10) 

This is the first order correctability condition. This means that with a suitable choice of ct and cy, a first order 
propagator can always be corrected to second order. Hence, the trace of any first order propagator is always second 
order. For example, the trace Tr(e^"^e'^^) is second order despite its appearance. 

With the first order correctability condition satisfied, the remaining commutators in H2.9|l are either [T, [T, V]] or 
[V, [T, V]], and can again be corrected by adding to C the term eCi = ecTv[T, V]. Thus with 

C = Co + eCi = ctT + cvV + £CTy[T, V] (2.11) 

such that [ct — cv) = ~^tv, we have 

Ha =T + V + e\eTTv[T, [T, V]] + evTv[V, [T, V]] ) 

+e^[Ci,T + V]+ e^erviCo, [T, V]] + ^e^[Co, [Co, T + V]] + 0{e^), 

= T + V + e'^icTTV - CTV + ^CTeTv)[T, [T, V]] 

+s^{evTV - CTV + \cvCTvW. [T, V\] + 0(2^) (2.12) 

If we now choose ctv = cttv + \ctCtv to eliminate the error term [T, [T, 1/]], then the error term [y, [T, 1/]] can 
vanish only if 

Cttv = cvTV -^{ctvY ■ (2-13) 

This is the general second order correctability requirement for correcting any first order propagator beyond second 
order. The major result of this work is to show that this condition cannot be satisfied for product decomposition of 
the form 12. 2() with only positive ti coefficients. 

III. DETERMINING THE ERROR COEFFICIENTS 

To check whether the correctability requirement H2.13|l can ever be satisfied by an approximate propagator of the 
product form (|2.2|) . we need to determine ctv ^ cttv and cvtv in terms of From the assumed equality 



(3.1) 



4 



with Ha given by (|2.2(l . we can expand both sides and compare terms order by order in powers of e. The left hand 
side of H3.1|l can be expanded as 



getiTgS.iVget.Tge..V ...get«igS.„V ^ I + , ^ y j^^ T + 6 \ } ] V, ] V + ■ ■ ■ , (3.2) 

and the right hand side as 

gSffA ^i + £{T + V)+ e^erviT, V] + e^eTTv[T, [T, V]] + e^evTv[V, [T, V]] 
+ ls^{T + V)^ + ^e^erv {{T + V)[T, V] + [T, V]{T + V)} 

+ ^e^T + Vf + ■■■ (3.3) 
Matching the first order terms in e gives the primary constraints 

N N 

^ti = l and ^v, = l. (3.4) 

i=l 1=1 

To determine the error coefficients, we "tag" a particular operator in (|3.3|l whose coefficient contains ctv, sttv or 
evTV and match the same operator's coefficients in the expansion of (|3.2|l . For example, in the terms of (I3.3|l . the 
coefficient of the operator TV is (^ + ctv) Equating this to the coefficients of TV from (|3.2|) gives 

1 ^ 

-+eTV = ^SiVi. (3.5) 

i=l 

where we have introduced the variable 

i 

= (3.6) 

This way of computing TV from H3.2|l corresponds to first picking out a V operator from among all the Vi terms, then 
combine all the ti terms to its left in the exponential to generate a T operator. Alternatively, the same coefficient can 
also be expressed as 

1 ^ 

+ evT ^^Uui. (3.7) 

■i=i 

where 

N 



This way of computing TV corresponds to first picking out a T operator from among all the U terms, then combine 
all the Vi terms to its right in the exponential to generate a V operator. To demonstrate how these variables are to 
be used, we can directly prove the equality of H3.5|l and (|3.7|l . First, note that sn = 1 and ui = 1. Second, since 
t'i — Si — Si-i, at z = 1 we must consistently set sq = 0. Similarly, since Vi = ui — Mi+i, we must set uat+i = 0. 
Therefore we have 

N N N N 

SiVi = ^ Si{ui - Ui+i) = '^{si - Si^i)ui = ^ (3.9) 

z— 1 i— 1 z— 1 z— 1 

The determination of error coefficients is simplified if we pick operators whose expansion coefficients are easy to 
calculate. Matching the coefficients of operators TTV and TVV (note, not the operator VTV) yields 

11 1 ^ 1 ^ 

^ + ^eTV + eTTV = ^'^sfv^ = -"^{s^, - s^,_l)u^, (3.10) 

1=1 1=1 

11 1 ^ 

- + -eTV-eTVT =-'^Uu1. (3.11) 

2 — 1 
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IV. PROVING THE MAIN RESULT 

Using the expression for ervT from H3.11|) . the correctability requirement H2.13|l reads 

1 ^ 

-Y^Uul^a, (4.1) 



with 

<^^\{\ + c^Tvf + Y^-eTTV (4.2) 

and ctv, sttv given by (|3.7|) . (|3.1UI) respectively. In Suzuki's prooS, he recognizes that in terms of the variable 
y/tiui, 14.1(1 is a hypersphere and (|3.7() . (|3.10|l are hyperplanes. His proof is based on a geometric demonstration that 
his hyperplane cannot intersect his hypersphere. While this geometric language is very appealing, it is cumbersome 
when dealing with more than one hyperplane. We will use a different strategy. 

If ti are all positive, then the LHS of (|4.1|l is a positive-definite quadratic form in Ui. There would be no real 
solutions for Ui if the minimum of the quadratic form is greater than a. Our strategy is therefore to minimize the 
quadratic form subject to constraints 1)3. 7(1 and 1(3.10(1 

N 

X^i^uz = b, (4.3) 

i=l 

N 

'^ti{si + Si^i)ui = c, (4.4) 



with 6 = i + evT, c = ^ + Ctv + '^sttv, and show that the resulting minimum is always greater than a. (The primary 
constraints 1(3.4(1 are just sat = 1 and ui = 1.) 

For constrained minimization, one can use the method of Lagrange multiplier. Minimizing 



1 ^ / w \ / N 

2=1 \i=l / \i=l 



ti{si + Si^i)ui - c] (4.5) 



gives 

w» = Ai + A2(s, + s,_i). (4.6) 
Substituting this back to satisfy constraints 1(4.3(1 and 1(4.4(1 determines Ai and A2: 

Ai + A2 = 6, (4.7) 

Ai + A2 + 5A2 = c. (4.8) 
The only non-trivial evaluation is X^il^i ^ii^i + Si-i)^ = 1 + g, where 

N 

9 = Y.'^sjs,-i-s,sj_i). (4.9) 

i=l 

The minimum of the quadratic form is therefore 

N 

2 

1=1 

[{Xi + + gXl] 



1 ^ 



i=l 

\[b' + -{c-bn (4.10) 
^ 9 
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To minimize F, one must maximize g. Solving dg/dsi = gives s; = (si+i + Si_i)/2, which means that Si is hnear 
in i. The normahzation sjq = \ fixes Si = i/N, giving 

ffmax = ^(l-^). (4.11) 

This is indeed a maximum since one can directly verify that d^g/dsf = — 2(si+i — Si_i) < 0. Hence, at any finite TV, 

F > + 3(c - bf] = + CTvf + ^(2 eTTy ^ = « + 6e|^^. (4.12) 

Thus the minimum of the quadratic form is always higher than the value required by the correctability condition. 
Hence, no real solutions for m are possible if U are all positive. 

We note that the above proof is independent of erv- For ctv = 0, the correctability condition is just cttv = evTV- 
Hence for symmetric decompositions with positive i^'s, where ctv = is automatic, we have as a corollary that bttv 
can never equal to cvtv- 



V. CORRECTABLE FORWARD PROPAGATORS AND THEIR CORRECTORS 



The last section is the main result of this work. Here, we show how the correctability criterion can be applied 
systematically to deduce forward correctable second order propagators and their minimal correctors. 

The proof of non-correctability is limited to the conventional product form (|2.2|) . which factorizes the propagator 
only in terms of operators T and V. As shown in the last section, symmetrically decomposed positive-time-step 
propagators cannot be corrected beyond second order because cttv cannot be made equal to evTV- For example, 
the second order propagator 

exp(ier)exp(£T/)exp(ier) (5.1) 

has ti = t2 = 1/2, vi — ui — 1, si = 1/2 and erv = 0. From (|3.10|) and H3.11|l . we can determine indeed that the 
two error coefficients are not equal: 

1 /1\^ 1 1 



A simple way to force them equal is to directly incorporate either operator [T, [T, V]] or [V, [T, V]] in the factorization 
process. Since [V, [T,V]] = l^i "(''uOP j^^* another potential function, Suzuki24 suggested that 

one should keep the operator [V, [T, V]]. If now we add ^£^[V, [T, V]] to eV in H5.1|l . we can change the coefficient 
evTV from —1/12 to —1/24, matching that of cttv- The result is still only a second order propagator 

PTi = exp QeT^ exp (^eV + ^e^'iV, [T, V]]^ exp QeT^ , (5.3) 

but now has a fourth order trace. This propagator was first obtained by Takahashi and Imad a^^i^^ by directly 
computing the trace. It is a remarkable find given how little they had to work with. This derivation explains, without 
doing any trace calculation, why the propagator worked. 

The alternative of keeping [T, [T, V]] would require adding -^£^[T, [T,V]] to make 

cttv equal to ey^y's value 

of —1/12. This operator is too complicated for practical use, but in the case of the harmonic oscillator, it can be 
combined with the kinetic energy operator: 

P2B = exp QeT - ^s' [T, [T, V]]^ exp (sV) exp QeT - ^e' [T, [T, V]]^ . (5.4) 

This can also be written in the form of 

P2B = exp (^ev\ exp (eT - ^^e^[T, [T, V]]\ exp {^ev\ . (5.5) 
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In this case exp{^eV) exp{sT) exp{^sV) has Cttv = 1/12 and eyrv = 1/24 and propagator p2B corresponds to 
changing cttv^^ value to match that of evTV- The Takahashi-Imada propagator (|5.3|) can also be written as 



p'ti = exp i^-eV + -s'[V, [T, V]] j exp (sT) exp l^-eV + -s'[V, [T, V]] j , (5.6) 

corresponding to changing e\/Tv's value to match that of cttv- These are the four fundamental correctable second 

order propagators with a fourth order trace. 

For the computation of the trace, it is unnecessary to know the corrector explicitly. In other cases, such as symplectic 

corrector algorithms, one may wish to apply the corrector occasionally to see the working of the corrected fourth order 

propagator p. We will give a detailed derivation of correctors for propagators (|5.3() - (|5.6() . cumulating in a set of four 

minimal correctors. These minimal correctors with analytical coefficients have not been previously described in the 
Iiteraturei8-i9-20.2i.22.23_ 

For the Takahashi-Imada propagator, we have cttv = &vtv = ^2 with 62 — —1/24. From l|2.7|l . we see that a 
possible corrector is C = e2e[T, V]. This can be constructed in a straightforward manner as suggested by Wisdom et 
alJ^. Since 

B{vi,ti) = exp{eviV) exp{etiT) exp{—eviV) exp{—etiT) 

= exp (^-v,t,e^[T, V] ~ ^t^ie^iT, [T, V]] - ^hvle^[V, [T, V]] + 0{e^)^ , (5.7) 

by setting viti — (1/48), the following product is a workable corrector 

B{vuti)B{-v,,-ti) ^ exp (^-^s^[T,V] + 0{s^)^ . (5.8) 

Note that it is important to have the operator V before T to generate a negative 62 coefficient. However, without 
fully determining both vi and ii, this corrector clearly under- utilizes B{vi,ti). It requires eight operators, which is 
far from optimal. We will show below that four is sufficient. 

Let H = T+V and G = [T, V]. Since Ha = H + C2e'^[H, G], we can see from (|T7j) that adding a term cqH to G, wiU 
not affect the corrector term e[G,T + V], but such a term will generate unwanted third order terms CQe2£^[H, [H, G]] 
from e[G,HA] and ^coe2e^[H, [G,H]] from ^s[G, [G,Ha]]- To cancel them, we must add another term C2e^[H,G] to 
the corrector such that C2 = ^0962. Thus the corrector can have the more general form 

exp(£C) = exp (^coeH + e2e^G + ^coe2e^[H, G]j + 0{e^), (5.9) 

= exp(coei?)exp(e2e^G) + 0(£''), (5.10) 

where the second line follows from the fundamental Baker-Campbell-Hausdorff formula, exp{A) exp{B) = exp(y4 -|- 
B + {1/2)[A, B] ■ ■ •). To exploit the use of the free parameter cq, we can approximate exp(coe-ff) by 

exp(ey F) exp(ecoT) exp(ey V") 

= exp (^cosH + ^cle'[T, [T, V]] + ^cle'[V, [T, V]]^ + 0(£^), (5.11) 

and the term exp(e2e^C?) by B{vi,ti). We can now choose co,wi,ii such that viti = 1/24 and the third order terms 
in (|5.11|l exactly cancel the third order terms in (|5.7|l : ^tfvi — j^Cq, ^tivf — ^Cq. This gives co = 1/(2 • 3^^^), 
fi = 1/(4V3) and ti = 1/(2^/3). The result is a corrector with six operators: 

S = exp(eyF) exp(£cor) exp(e(y + vi)V) exp(eiiT) exp(-ewiV) exp(-eiiT). (5.12) 

Since this corrector has made good use of all the parameters, it is surprising that one can find a even shorter corrector. 
Instead of B(vi, ti), consider just 

exp{edoV) exp{edoT) 

= exp UeH - idg[T, V] + ^dy[T, [T, V]] - ^dy[V, [T, V]]\ + Die'). (5.13) 
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The corrector 

Sti = exp{e^V)exp{ecoT)exp{e{^ + do)V)exp{edoT) (5.14) 
= exp(^(co + do)eH + {-^dl)e'G + ^i-^dDico + d„)[H, G] 
+ ^(4 + ^dDs^T, [T, V]] + 2^(cg + Ml)e'[V, [T, V]]^ + 0{e') 

will have the correct value for 62 if we take o?o/2 = 1/24, fixing do = l/(2\/3). The corrector will also be of the 
form 1)5. 9(1 after both commutators have been eliminated by setting Cq = — 4(iQ, giving cq = — 1/(2^/'^\/3)- This is the 
minimal corrector for the Takahashi-Imada propagator. 

The corrector of the form l|5.9|l is completely determined by a single number 62. Its sign dictates the order of the 
T and V operators and its value fixes their coefficients. For the alternative propagator H5.4I) with 62 = —1/12, 
its corrector is of the same form as l|5.14f) . but now with c?o = l/\/6 and co = — 2^/^/\/3. 

For positive values of €2 , the corrector is of the form 

S = exp(eyr)exp(ecoF)exp(e(y + do)r)exp(edoV") (5.15) 
= exp(^(co + d^)eH + {\dl)e^G + \{\dl){ca + da)[H, G] 

^Ya^cI + ^dl)e^[T, [T, V]] - ^{cl + Ml)e'[V, [T, V]]^ + 0{e^). 

Propagator p2B is dual to the TI propagator with 62 — 1/24. Its corrector is of the form (|5.15|) but with same 
coefficients do = l/(2-\/3) and cq ~ —1/(2^/'^ a/3)- The p^j propagator H5.6() with 62 = 1/12 is dual to P2g. Its 
corrector is of the form (|5.15|) with do = 1/V6 and co — — 2^/^/\/3. These compact correctors are fitting companions 
to their equally compact propagators. 

VI. CONCLUSIONS 

In this work, we proved a fundamental result on the correctability of forward time step propagators. We show that 
if p = e'^(-^+^^ were to be approximated by the product form (|2.2|l . then no product form with positive coefficients {ti\ 
is correctable beyond second order. Whereas a conventional higher order propagator requires its error terms to vanish, 
a correctable propagator only require its error terms to satisfy the correctability condition. The latter requirement 
seemed far less stringent. A surprising element of this work is that, this is not the case. For symmetric decomposition 
with positive {ti}, the two second order error coefficients cannot both vanish because, they can never be equal! The 
correctability requirement itself is stringent enough. This proof of non-correctability generalizes the previous work of 
Shengl and Suzuki^. 

From knowing correctability requirement, we derived systematically the four forward correctable second order 
propagators and their minimal correctors. These minimal correctors follow from a more general form H5.10|l of the 
corrector with free parameters. Much of the existing literature on symplectic corrector is rather opaque, concerned 
only with how to satisfy "order conditions" numericallyS2i2^. This work suggests that a more analytical approach is 
possible. 

The Takahashi-Imada type of propagators considered here are unique in that they are the only known second-order, 
forward-time-stcp propagators with a fourth order trace. If one is willing to evaluate the potential at least twice, 
then with the inclusion of \V, \T,V]], one can make both error coefficients cttv and evTV vanish^'--. The result is 
a whole family of positive time step fourth order propagator32L22iS2i^, with a fourth order trace. While this class 
of forward decomposition algorithms is indispensable for solving time-irreversible problemsi^iiLiSiiiii, they are less 
interesting from the point of view of calculating the trace. For correctable propagators, their key attraction is that 
one can obtain a higher order trace without using a higher order propagator. Methods and results of this work can 
be used to study ways of correcting these fourth order propagators to higher orders. 
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